Gap probability of a one-dimensional gas model 



Y Chen 1 ^ and S. M. Manning 2 ^ 
1 Department of Mathematics, Imperial College 
180 Queen's Gate,London SW7 2BZ UK 
2 Department of Theoretical Physics, Oxford University 
1 Keble Road,OXl 3NP UK 
February 1, 2008 

Abstract 

We investigate the gap formation probability of the effective one dimensional 
gas model recently proposed for the energy level statistics for disordered solids at 
the mobility edge. It is found that in order to get the correct form for the gap 
probability of this model, the thermodynamic limit must be taken very carefully. 
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Introduction 

The eigenvalue distribution of large random matrices has a long history beginning with the 
work of Wigner Q on the energy level statistics of heavy nuclei. This topic has found recent 
applications in the theory of disordered metals, quantum chaos and string theory. The most 
studied and best understood of random matrix models are the Gaussian ensembles with unitary, 
orthogonal and symplectic symmetries ||]. Perhaps the most fundamental of all the statistical 
quantities of interest is the probability that an appropriately scaled interval of the spectrum is 
free of eigenvalues. This is a measure of the degree of correlation between the eigenvalues. There 
is an ansatz of Wigner which gives a simple formula that fits rather well the experimentally 
observed nearest levels spacing distribution § Q. The energy levels of disordered metals in 
the weak localization regime exhibit behavior which is rather well described by the correlation 
functions of the various Gaussian ensembles || (7|. Furthermore the associated level spacing 
distribution function compares favorably with the Wigner Ansatz, but upon increasing the 
strength of disorder the level spacing distribution evolves and deviates from the "standard 
behavior." Perhaps the most striking result is that at the mobility edge the level spacing 
distribution has a universal form when expressed in units of the average level spacings ||. A 
random matrix model has been proposed which accounts for this continuous deviation Q. In 
this model the usual unitary Gaussian ensemble is q— deformed in a similar way that to the 



q— deformation of the Laguerre ensembles [ 10 1 . The results of this ensemble give qualitative 



features which are in agreement with the energy level statistics obtained from numerical work on 
disordered solids [§]• Using the two-level correlation function obtained in 0] it was established 



1 Although this ansatz deviates only minutely from results of numerical experiments on random matri- 
ces, it is now known that the Wigner ansatz is not quite correct; the complete level spacing distribution 
functions are obtained from transcendental functions of various Painleve types; the type depends the part 
of the spectrum where one performs the scaling || [|, |j| 



in [11] that the number variance and the A3 statistics exhibit novel behavior. Very recently, 



using a combination of scaling arguments and diagrammatical analysis, it has been suggested by 
the authors of that the density of state correlation function in the bulk of spectrum at the 
mobility edge behaves as —\/\x — y| 2-7 instead of — l/|x — y\ 2 in the regime of weak disorder. 
Here x, y denotes the energy levels and 7 is related to the correlation length exponent v and 
the dimension (d) of the system: 7 = 1 — -K. 

A gas model — known as the effective plasma model in [14, [0J — with a power law pair repulsion 



between the particles which gives correlation function: — l/|x — y\ 2 7 , was proposed. Using a 



continuum approximation of Dyson [16], the authors of [14, [15|] conclude that there is a new 
universal law for the level spacing distribution function for the energy levels of a disordered 
metal at the mobility edge. This law states that the probability that the interval (— s, s) of the 
spectrum is free of levels behaves as E(—s,s) ~ e -constant s 2 1 ^ where the constant appearing 
here is independent of s. 

In the continuum limit, which is expected to be valid in the limit where N (the total number 
of particles) goes to infinity, the gas model may be described by a continuous charge density a, 
following the work of Dyson . In the mean field approximation this gas model has the energy 
functional, 

F[a] = ~ / dx [ dy ^^ + / dxa(x)u(x) , < 7 < 1 (1) 
2 J j J j \x-yp J j 

subject to 

J dxa{x) = N, 

where J is the support of the eigenvalues/fluid and u{x) is the confining potential that holds 
together the repelling particles. Without loss of generality the coefficient of l/|x — y\" ( in Eq. 



(1) is set to unity. This is denoted by Aj in [14, 15]. Although the continuum approximation 



is not mathematically rigorous, previous experience suggests that it is rather robust p2j. The 



density that minimizes Eq.(l) satisfies, 

J dy-^ = A — u(x), x € J, A = Chemical potential (2) 

supplemented by Jj dxa(x) = N. At equilibrium the minimum free energy is 

F[J] : = F chem + F int , (3) 

where -F c hem := AN/ 2 is chemical potential contribution to the free energy and 
-^int '■= 2 fjdxa(x)u(x), is the energy due to interaction between the charge distribution and 
the confining potential. Eq.(3) can be obtained simply by the use of the mean field equation 
that governs a and the normalization condition. Note that the separation into -F c hem and Fi nt 
in Eq.(3) is not unique, different but equivalent expression can be found by using Eqns. (1) and 
(2) and the normalization condition. To obtain the density we have to deal with the inversion 
of the following integral equation: 

/ _^Ld y = v>(x), xe J. (4) 
J j \x - yp 

where < 7 < 1, and J is a subset of the real line. Assuming i\) £ L 2 [J], we seek (p € £ 2 [«7]. In 
the case where J is a single interval, say J = (0, o), there is a unique inversion formula, 

B y d f a , t 1 '^ d ft , Ma^r 



<P( X ) = ^=T^Z dt ~ TEE ds Z -ttt", x G (0,a), (5) 

x ^dxJ x (t-x)^ 1 dt Jo (t-s)^ 

where 

r( 7 )cos(=) 



7 

7T 



r(l?)l 2 



Eq. (4) was investigated by Carleman [17], and later by Widom in the context of stable processes 



[18|. This equation is also a special case of an integral equation with hyper geometric kernel [19]. 



However, in the many interval case, 

J = Up =1 (o p , b p ), ai < bi < . . . < a n < b n , 
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no similar result has been obtained. This is an outstanding problem. 

The quantity of interest, the probability, E[J], that an interval J is free of eigenvalues is 

E[J] _ (n^i^^)e- 

where 



-W 



rij=l fju J c dxj J e 



N 1 



j=l l<j<k<N |X J Xfc ' 7 

J c is the complement of J and J c U J is the natural support of the particles. 
From the previous equation we see that minus the logarithm of E[J] is the change in free energy; 
the free energy where all particles reside in J c minus the free energy where all particles reside 
in J c U J, 

-In E[J] = 5F:=F[J C }-F[JUJ C }. 



The model considered in [14, 15 1, corresponds to the particles being distributed on the real line 
and and has the confining potential u(x) = x 2 . The logarithm of the probability that the interval 
(—a, a) is free of particles is formally 

- lnE[(-a, a)] = 5F := F[(-b, -a) U (a, b)} - F[(-B, B)] 

= \n [A(a, N) - A(0, N)} + J / dxa(x, gap)n(x) - \ f* a(x)u(x), (6) 

I I J(-b,-a)U(a,b) 2, J-B 

where b(> a > 0) is band edge of the spectrum with a gap at (—a, a) and B is the band edge 
where the spectrum is gapless. We note that with the normalization condition implemented, 
the respective edges B and b are functions of N, and will diverge as N — > oo. Furthermore, b 
is also a function of a. The deliberate use of B and b to distinguish the two situations is to 
indicate that although in the thermodynamic limit, i.e. N — > oo, B and b differs by a very small 
amount their difference is nevertheless important and can not be neglected, for in computing 
5F we are subtracting two very large numbers. Concerning the chemical potentials; A(a, N) 



is the chemical potential where there is a gap and A(0, N) is that without. The above remark 
also applies to the change in the chemical potentials. The density, cr(x,gap), indicates that the 
integral equation must be inverted in the interval (—6, —a) U (a, b) for which an explicit inversion 
formula is not known. 

Although the change in the chemical potential is expected to be small as N —> oo but due 
to the multiplicative factor N; this will give a significant contribution to the change in free 



energy. The authors of [14, 15] assume that the change in the chemical potential (denoted 



by /i in [14, [15|]), Sfj,, is of order s/£, (where their s is our a and £ is our B) and vanishes 
in the thermodynamic limit, because £ diverges as N — > oo, giving negligible contribution 
to 8F. However, further reflection suggests otherwise. From the symmetrical density given in 
[14, [15| , p{x) ~ (£ 2 — x 2 )^ 1 , a simple calculation using the normalization condition shows 

1 1+7 

that £ ~ N 2 +~>, thus (5-Fchem ~ NSfi ~ N 2 +~<s and cannot be neglected. Indeed the scaling 
variable S commonly used in random matrix theory ( not the s used in 14, contains the 
product of two quantities one of which — s — tends to zero and the other a function of N tends to 



infinity. Therefore, the appropriate scaling variable for the model studied in [14, 15] ought to be 



1 + 7 

S = N 2 +f s. In an example to be given below, where the problem of multi-interval inversion of 
the integral equation is circumvented, the meaning of the scaling variable S will be clarified in 
an exact computation within the mean field approximation. In fact we will see that the change 
in chemical potential is of order a, and is not of order a/B. 
Gap formation probability of a semi-infinite model 

As mentioned in the previous section in order to by-pass the problem of multi-interval inversion, 
we consider the model with u{x) = x, x > 0. We are required to compute two free energies; one 
with all the particles are confined in the interval (o, b), where b > a > 0, and the other with all 
the particles confined in (0,B). Here b is the band edge for a density which has a gap in (0, a), 
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and is supported only in (a, b). The band edge B is for a density that is supported over (0, B). 
Note that when the normalization condition is implemented b will be a function of N and a 
while B will be a function of N only. The probability that the interval (0, a) has no particle in 
the thermodynamic limit is therefore 

E[S] = exp [-SF] , (7) 

where 

SF := F[(a, b)] - F[(0, B)]. (8) 

The scaling variable S will become apparent later. To determine the free energy, we first invert 
the integral equation, 

[ b dy-^- = A-x, a<x<b, (9) 
J a \x-yy 

to find a. With a shift in the variables x — ► x + a, y ^ y + a we can now use Eq. (5) to find, 
for a < x < b, 



a[x) — 



(x _ a) V r( 7 ) 



(10) 



2 7 / (ft-x) - ^ 7 
Note that since < 7 < 1 the density will have integrable divergences at a and b. In order to 
have a configuration that has a lower free energy we choose 

A = {b-a)(^)+a (11) 

to eliminate the divergence at b. This is in agreement with the general theory of integral equations 
with weakly singular kernel. The normalization condition, dxa(x) = N, produces, 

C o S nn.) r 2 (lp) 

n = — n 2 ' \ J (b- a y +1 . (12) 

2vT7 2 r(7) v ; v ; 

From Eqns. (11) and (12), we find 

'ZE2) p 2 ^ 1 +" 
Cjl ' ~ 7 '~ 2vr7 2 r( 7 ) 



and 



7 + 1\ IN 



i 



1 

Eq. (13) shows that I? = (_/V/C 7 )t +1 . The first term in the r.h.s. of Eq. (14) is identified as 
A(N,0), the gapless chemical potential. The change in the chemical potential is therefore a and 
not a/B as claimed in |l4| , [l5| . With a given by Eq. (10) and A(N,a) given by Eq. (14), we 
compute the free energy in the presence of the gap, using elementary integration formulae, 

Fmt ~ — + 2vr 7 r( 7 + 3 ) [a;) ' (15) 

7 + 1 / N \ ^ Na 

F ^-— N {c;) < 16 > 

Therefore the free energy for which all particles are confined in (a, b) is 

F(a,b)=Na + f 7 {N), (17) 

where 



T+2 

7 + 1 / N ' 



We identify f-y(N) with F(0,B), the free energy for which all particles are confined in (0, B) 
and the change in free energy is exactly, 

5F(0, a) = F(a, b) — F(0, B) = Na. (19) 

The limit which produces the scaling variable is defined as a — ► and N — > oo such that S := Na 
is finite. Thus the probability that the interval (0, a) is free of particles is 

E[(0,a)] = e' s . (20) 

To further illustrate the meaning of 5, we consider another model with u(x) := x 2 , x > 0. Using 
Eq.(5), and after some computation, we find for a < x < b, 



cos (El) 1 r 
a ( x ) = ~ 7ZI7I~TiT ~ TT=7 C+[4a- 12b - 4(o + 6)7] (b - x) + 8(6 - x) 



4^7(7 + 1) [( a ._ a )( iC _6)]T : 



(21) 



where 



C := -a A - 2ab + 3b z - 4A 7 ( 7 + 1) + 46 z 7 + (a + 6) 2 7 2 . 



(22) 



Just as in the previous example we choose 



A 



(a + 6) 2 7 + 3b 2 -a 2 - 2ab 
4 7 



(23) 



so that C = 0, to eliminate the integrable singularity of a at 6. Thus 



a ix) 



cos g) 



, 1 + 7 
X) 2 



x— [( 7 + 3)b - (1 - 7 )a - 2(6 - x)] , a < x < b. 



(24) 



^7(7 + 1) ( x -a)V 

We note that the factor [. . .] in Eq.(24) for a < x < b is greater than (a+6)( 7 + l) and is therefore 
positive. Application of the normalization condition leads to a transcendental equation for b, 



N 



D^ 2 (jg 

2r( 7 + i) 



(b-a) 



7+2 



1 , 7(^ + «) , 7 + 1 



b-a 7 + 2 



7T 7 ( 7 + 1) 



(25) 



Although the explicit solution of the transcendental equation governing b is not known we may 
nevertheless determine 6 as a power series in a, which in turn can be used to calculate Fi nt 
and -Fchem, also in a power series in a. After some straightforward but lengthy calculations not 
reproduced here, we find 



5F ~ aiv'S 1 + 



In this case the appropriate scaling variable is 



(26) 



3+7 

S ■= o# 2 +t, 



(27) 



and 



E[(a,b)\ ~ e" 



(28) 



We see that although the scaling variable S is the same in both cases, the thermodynamic limit 
is arrived in a distinct manner for distinct potentials. Following the procedure described in 
we will find — lni?[s] ~ s for our potentials. To see how this comes about, we compute the 



9 



change in the free energy according to the procedure given in [15| for the semi- infinite models 



studied above. In order to determine the probability that there are no particles in (0, s), we 



require the change in free energy, which according to [15|, is 



F . - Fo . -i r * r ix .s^gi + \ r dx r (M) 

2 Js Js \x — x'p 2 Jo Jo \x — x'p 

where, 5p s (x), the change in the density due to the presence of the gap is related to the density 
in the absence of the gap, po(x), via the following integral equation: 

r dx .^^=r dx . T M*L. ( 30) 

Js \x-x'\~' Jo \x-xf\i K ' 

Eq. (30) which governs 5p s and which provides a determination of F s — Fq is equivalent to the 



statement that the change in the chemical potential vanishes in the thermodynamic limit [15]. 
However, as noted from the examples given above that the change of chemical potential is of 
order a. We now deduce the s dependence of the gap formation probability from Eqns. (29) and 
(30). With the change of variables used in fll5|| : z = s/x, t = s/x', and u g (z) = \z\~ 2+1 5p s {s / z), 
Eq. (30) becomes, 

I <B j^ = 4 T r (31) 



I 7 s ~ •'i t~\z - t[y 

To obtain this we have used the fact the integral in the r.h.s. of Eq. (30) is dominated by po{x) 

i — f 

near the origin: po(x) ~ l/x^~ . Note that the t integration in the r.h.s. of Eq.(31) converges. 

1 — 7 

From the Eq. (31) we see that the combination u g (t)s^~ is independent of s. Using this, and 
the variables z, t and u g (t), we find from a simple dimensional argument, 

F s — Fq = constant s, 

where the constant is independent of s. This suggests that the change in the free energy is 
independent of the confining potential u(x) and universal. The statement of universality requires 
some qualification. If the confining potential u(x) is far weaker then those considered above, the 

10 



singularity of density near the origin, (which is determined uniquely by the confining potential 
via the integral equation,) will be modified and we may expect the final form of ln.E[,s] to be 
modified accordingly. 

It is interesting to note that for sufficiently strong confining potentials the singularity of the 

1-7 

density is universal, i.e., po(x) ~ 1/x 2 . This phenomena is also observed in the standard 
random matrix ensembles delineating the behaviour of the density for strong and weak confining 
potentials (no longer universal) (l(], |9|, [ll], p^| . Thus in the situations where the singularity 
structure of the density is universal a universal form for — In E[s] is found. 
It appeared by adopting the two distinct approaches; one of which computes the change in 
the chemical potential explicitly by keeping N finite but large while the other taking the limit 
N — > oo right off the start, we are led the same law for the gap probability, but with distinct 
scaling variable S = N v a as opposed to s. It may be of interest to understand better the 
thermodynamic limits and clarify the meaning of s and S. 

We are at present developing tools in an attempt to obtain a multi-interval inversion formula. 
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